Andreev Level Spectroscopy and Josephson Current Switching in a 3-Terminal 

Josephson Junction 
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We calculate the electrical currents through a superconductor - insulator - superconductor junc- 
tion which is also weakly coupled to a normal metal side probe. The voltage V applied to the 
normal metal terminal controls the occupation of Andreev energy levels En, and therefore controls 
the Josephson current flowing through these levels. Whenever the probe voltage crosses an Andreev 
level, the Josephson current changes abruptly by an amount equal to the current flowing through the 
Andreev level. The differential conductance along the normal metal terminal permits spectroscopy 
of the Andreev levels. In a short junction {L ^ ^o), the critical current switches abruptly from 
the Ambegaokar-Baratoff value to zero when the probe voltage is approximately equal to the su- 
perconducting energy gap {\eV\ ~ A). The magnitude of the Josephson current switching in a long 
junction (L 2> Co) , and the range of probe voltages over which the Josephson current differs from 
its equilibrium value, are much smaller than for three-terminal ballistic superconductor - normal 
metal - superconductor junctions. 

PACS numbers: 74.80Fp, 74.50-Fr, 73.20.Dx 

submitted to Physical Review B-1 

I. INTRODUCTION elastic processes is extremely small, even a normal metal 

probe only weakly in contact with the Josephson junction 
will inject quasi-particles into the Andreev levels much 
faster than the superconducting electrodes. As long as 
the normal metal probe is only weakly in contact with the 
Josephson junction through a tunnel barrier, the main ef- 
fect of the probe is therefore to fix the occupation factor 
of the Andreev levels, leaving the wavefunctions of the 
Andreev levels essentially unchanged from the isolated 
Josephson junction.Elllj'Ej The normal metal probe there- 
fore forces the effective Fermi level of the bound Andreev 
states towards the Fermi level of the probe, rather than 
the Fermi level of the superconductors. 



The Andreev energy levelso in superconductor - normal 
metal - superconductor (SNS) or superconductor - insu- 
lator - superconductor (SIS) junctions, through whji 
a large fraction of the the Josephson current flowsa" 
are only weakly held in thermodynamic equilibrium with 
the two superconducting contacts of a Josephson junc- 
tion. When the quasi-particle energy is inside the su- 
perconducting gap, quasi-particles cannot transmit into 
the superconductor, and also cannot be iiijected into the 
Andreev levels from the superconductom. The coupled 
electron-like and hole-like quasi-particles which form the 
Andreev levels orbit in continuous periodic motion in- 
side the normal (or insulating) region of the SNS or SIS 
junction. Quasi-particles in the Andreev levels are then 
essentially thermodynamically isolated from the super- 
conductors, yet carry a large fraction of the supercur- 
rent. Only inelastic scattering inside the superconductor 
forces the occupation factor for the Andreev energy levels 
towards the equilibrium Fermi occupation factor of the 
superconducting contacts. 

An additional normal metal contact coupled to the 
Josephson junction can directly inject quasi-particles into 
the Andreev energy levels through elastic scattering pro- 
cesses, and therefore can directly control the occupation 
of the bound levels. The additional normal metal probe 
coupled to the Josephson junction, shown in Fig. |^, mod- 
els either a scanning tunneling microscope tip or the_ffa±e 
electrode of a three-terminal Josephson junctionE3il3. 
Since the rate at which the superconducting contacts in- 
ject quasi-particles into the Andreev levels through in- 
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FIG. 1. Josephson junction coupled to a normal metal side 
probe. The probe is biased at a voltage V with respect 
to the two superconductors. This normal metal probe con- 
trols the occupation factor of the bound Andreev levels, since 
quasi-particles having energies inside the superconducting gap 
cannot transmit into either superconductor. 

Controlling the Andreev bound state occupation 
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through a normal terminal leads both to an abrupt 
switching of the Josephson current I{V), and a peak in 
the differential conductance dI]\[{V)/dV along the nor- 
mal terminal, whenever the pipbe voltage V is equal to 
the energy of an Andreev level.ll3 As a bound level is pop- 
ulated or depopulated by the probe voltage, the Joseph- 
son current changes by an amount equal to the current 
carried by the bound level. As the Andreev level is being 
filled by the probe, the differential conductance diN /dV 
along the normal metal lead also has a peak. The den- 
sity of Andreev levels in the SNS junction can therefore 
be detected by measuring the differential conductance 
diN/dV along the normal metal lead. Together with a 
phase - biasing network of large area Josephson junc- 
tions, this differentiaL conductance spectroscopy can be 
used to directly mapll3 the energy - phase relation En{(j)) 
and current - phase relation /„(</') of the Andreev energy 
levels in different types of Josephson junctions. 

One can obtain a more quantitative understanding of 
the Josephson current switching and Andreev level spec- 
troscopy in any Josephson junction in the limit where 
the normal metal probe is weakly coupled to the junc- 
tion. The total Josephson current /(</), V) is a sum of 
the currents flowing through discrete energies inside the 
superconducting energy gap Id and through the contin- 
uum of energy levels outside the gap Ic as /(</), V) = 
/d(</), V)+Ic{4>)- Here (/) = (5^2 — 0i is the superconducting 
phase difference. The contribution Ic to the Josephson 
current by scattering states outside the superconduct- 
ing energy gap is essentially unchanged by the probe 
voltage V, since the superconducting contacts can eas- 
ily inject quasi-particles into Andreev resonances in the 
energy continuum. Scattering states outside the energy 
gap therefore remain in equilibrium with the supercon- 
ducting contacts. The portion of the Josephson cprrent 
Id flowing through the Andreev levels, however, igl3 

V) = EnUn W f {En W " eV) + 

J+(</.)/(i?+(0)-ey)} • ^ ' 

In Eq. (1) the /±(<;i) = {2e/h){dE^{(j))/d(j)) are the cur- 
rents carried by 'forward' and 'reverse' Andreev levels 
E^ before adding the side probe. The probe voltage 
eV appears inside the Fermi factors / in Eq. (0) as an 
effective electrochemical potential for the Andreev lev- 
els. Therefore, Eq. (|l|) implies that the contribution of 
each Andreev level to the total Josephson current can be 
switched on or off by varying the probe voltage V. 

The tunneling current through the normal metal probe 
measures the local density of quasi-particle states in the 
Josephson junction. It is well known from tunneling spec- 
troscopy of normal metals that the tunneling current is 
propartionai to the local density of states at the sur- 
facella, and this also holds true for superconducting tun- 
nel junctions. The tunneling spectroscopy of Andreev 
levels in a Josephson junction, using tunneling current 
from the normal probe then, corresponds to 



dV h l^^{rl + ieV-E^^)^) ' ^ > 

when \eV\ < A. Here r„ is the width of Andreev level n, 
which is profiortional to the coupling constant e. Al- 
though Ref.L3 derived Eqs. only for a ballistic 
SNS junction, they should describe any type of Joseph- 
son junction. The numerical simulations we present in 
the following sections follow from the scattering theory 
in Appendix A, and can be understood using Eqs. 

In this paper we consider the Josephson current switch- 
ing I{V) and differential conductance spectroscopy of 
the Andreev levels dlN{V)/dV as we vary the voltage V 
along the normal terminal in a three-terminal SIS junc- 
tion. The details of the Josephson current switching and 
spectroscopy of the Andreev states in a three-terminal 
SIS junction differ considerably from the ballistic SNS 
junctional. In a short SIS junction, having (L <C (,o), 
the presence of an insulator forces the Andreev energy 
levels to the edge of the superconducting energy gap. 
The Josephson current in a short SIS junction therefore 
switches to zero when the voltage on the normal terminal 
is approximately equal to the energy gap, i.e. eV ~ A, as 
we discuss in section 3. Section 4 shows that the termi- 
nal I-V characteristics of SIS Josephson junctions longer 
than the healing length (L 3> Co) a-re more complex than 
those of short SIS junctions. The size of the nonequilib- 
rium Josephson current, the regularity of its switching 
behavior, and the voltage range over which the terminal 
currents are constant, are sensitive both to the barrier 
transmission T and the position of the tunnel barrier in 
the junction. We also discuss the special limiting ease 
where the long SIS junction has inversion symmetryc3. 



II. SHORT JOSEPHSON JUNCTION 



In a short SNS junction {L ^ ^o), the Josephson cur- 
rent flowing into either superconductor switches on or 
off as we vary the bias voltage on the side probe. It is 
well known both that short Josephson junctions contain 
only two Andreev levels, and that all the Josephson cur- 
rent flows through these levels {Ic = 0). Depopulating 
(or populating) both levels therefore forces the Josephson 
current to zero. Consider the SIS junction having trans- 
mission probabihty T = 2.5% in Fig. |. The two Andreev 
levels En{(j)) for the SIS junction are shown in Fig. ||(a). 
The horizontal lines in Fig. ^(a) correspond to bias volt- 
ages near the energy gap, namely eV = ±0.995A. For 
positive bias voltages, injection of a quasi-particle from 
the normal metal terminal fills any Andreev level hav- 
ing energy less than En < eV. When only the lowest 
Andreev level satisfies En < eV, only the lowest en- 
ergy level contains a quasi-particle. Consequently, when 
eV < VRA,U the Josephson current through the SIS 
junction is nearly the same as without the probe, i.e. we 
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recover the standard Ambegaokar-BaratofF resulti'lli~ii, 
as shown in Fig. ||(b). 
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FIG. 2. (a) Andreev levels E^{(l>) (solid) and bias voltage 
eV = ±0.995A (dashed) in a short SIS junction, (b) The 
probe voltage crossing an Andreev level forces the Josephson 
current to zero, except for a small leakage current. 



4> = 0.3(27r) (b)-(c). With decreasing transmission T, 
an energy gap opens and forces the Andreev levels to 
the superconducting gap edge in Fig. ||(a). Lowering the 
transmission coefficient in SIS junctions, therefore forces 
the switching voltage and differential conductance peak 
towards the superconducting gap at eV = A, as shown 
in Fig. ^(b)-(c), respectively. 
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By varying the gate voltage, we can switch the Joseph- 
son current in a short SNS junction on or off. At a 
fixed phase difference, the two Andreev levels carry equal 
amounts of current but in opposite directions. There- 
fore, when both Andreev levels are filled (or empty) the 
Josephson current is nearly zero as shown in Fig. 0(b). 
We also conclude that the Josephson current-phase re- 
lation is nearly the same whether the side probe has a 
negative or positive bias voltage. The small difference 
between the two Josephson currents for eV = ±0.995A 
shown in Fig. ||(b) is due to the small leakage current 
from the gate. The leakage current is small because the 
coupling strength e — 0.1% in Fig. |^. The small leakage 
current implies that, as we vary the gate voltage between 
VRA < eV < A, the Josephson current switches from 
the Ambegaokar-Baratoff value to approximately zero. 

Reducing the insulator transmission T forces the 
switching voltage and differential conductance peak to- 
wards the energy gap, i.e. eV — A, as shown in 
Fig. ^. In Fig. ^ the transmission probability decreases 
from T — 100%, 36%, 12%, the probe couphng is weak 
(e ~ 0.1% in (a)-(b) and e = 5% in (c)), and we fix 
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FIG. 3. (a) Andreev levels E^((j>), (b) Josephson current 
I{(j>,V), and (c) differential conductance diM /dV {(j},V) along 
the normal metal probe. Decreasing transmission probability 
(T = 100%, 36%, 12%) forces both the switching voltage in 
(b) and peak in the differential conductance in (c) towards 
the gap edge. 
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FIG. 4. A long (L > ,^o) SIS junction which has inversion symmetry, (a) Andreev levels, (b)-(c) Josephson currents, and (d) 
differential conductance along the normal probe. With decreasing barrier transmission, the nonequilibrium Josephson current 
exceeds the equilibrium current {V = 0) over a narrow voltage range. Ballistic junctions (T = 100%) have the largest currents 
both in and out of equilibrium. 



We can understand why the Josephson current in 
Fig. ^(b) is constant before switching off when the bias 
voltage crosses an Andreev level at £'„ = eV by con- 
sidering the Andreev level structure in Fig. ^(a). For a 
positive bias voltages and fixed phase difference 0, the 
lowest Andreev level (En < 0) carries the total current 
until the upper level crosses £'„ = eV. The zero temper- 
ature Josephson current is therefore unchanged for volt- 
ages smaller than \eV\ < |i?„|. Decreasing T reduces the 
Josephson current in Fig. ^(b) , since the lowest Andreev 
level carries a smaller current with smaller transmission 
T. The differential conductance (Un /dV along the probe 
also has a peak whenever the bias voltage crosses a new 
Andreev level, as shown by comparing Fig. ^(b)-(c). 

III. LONG JOSEPHSON JUNCTION 



Spectroscopy of the Andreev levels and the Josephson 
current switching as a function of the gate voltage V 
change significantly when the junction length L becomes 
comparable to the BCS healing length ^o- The number 
of levels is proportional to i/^Oj so more Andreev lev- 
els En{(j3) become bound in the pair potential well. In 
addition, the interference pattern between quasi-particle 
waves multiply reflected between the NS interface and the 
tunnel barrier depends on the junction length L/^o ^^nd 
the symmetry of the SIS junctioncH. The current In{4>) 
flowing through the Andreev levels also depends on the 
junction length and symmetry of the scattering poten- 
tial. To illustrate the variation of the Josephson current 
and probe current with these parameters, we study in 
this section a 'symmetric' SIS junction (where the impu- 
rity is in the middle of the normal metal region) and an 
'asymmetric' junction lacking this inversion symmetry. 
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A. Symmetric Junction 

It has long been known that the net equihbrium 
Josephson current in long SNS and SIS junctions is a 
small difference between much larger jBcisitiya, and neg- 
ative currents flowing in equilibriuinl3'CrE3i3i3. For a 
Josephson junction in equilibrium, adding these large 
counterflowing currents produces a net current propor- 
tional to the barrier transmission T. Ref.E3 quantified 
the magnitude of the larger electrical currents (which 
cancel in equilibrium) . In long, low transmissioa Joseph- 
son junction having inversion symmetry Ref-EZI showed, 
that these larger currents are proportional to ^/T. Ref.EZi 
also suggested that one could probe these 'giant' currents 
by inducing a non-equilibrium population of the Andreev 
levels. 

Fig. ^ shows the (a) Andreev levels, (b)-(c) Joseph- 
son currents, and (d) differential conductance along the 
normal probe in a long {L — 6.6^0)1 symmetric SIS junc- 
tion. The probe coupling in Fig. ^ is also weak, so that 
e = 0.1% in (a)-(c) and e = 5% in (d). The Andreev levels 
in Fig. |^(a) do not split at = 0, 27r, as the transmission 
decreases from T = 100% (fiUgd) to 36% (empty) or 12% 
(dotted), as one might expect. El The Andreev levels in the 
symmetrical junction do not split at (/) = 0, 27r because 
of a geometrical symmetry in the junction. Upon normal 
reflection from the tunnel barrier, a quasi-particle cannot 
tell whether it is on the left or right side of the barrier 
when = 0, 27r. Therefore, the energy levels at </> = 0, 2n 
are unaffected by the presence of a tunnel barrier. When 
the phase difference is not = 0, 27r this geometrical 
symmetry is broken, so the degenerate energy levels do 
split at any phase value other than = 0, 2n. This fail- 
ure of the energy levels to split at = 0, 27r in s-wave SIS 
junctions is exactly the same geometrical symmetry lead- 
ing to the 'midgap' energy levels in Josephson junctions 
formed from d-wave superconductor^3iHj. 

To obtain a formula for the Andreev levels in Fig. ^(a), 
we set a = L/2 in Eqs. (7)-(8) of Ref.l. The 'effective 
phase' a in the symmetric SIS junction then simplifies to 

sin(a) ~ a ~ 2\/T|sin(^)|, (3) 
leading to the Andreev levels 

Et ^ ;^f^(2^" - TT T 2^/T| sin(^)|) . (4) 

The geometrical symmetry also produces an additional 
resonant enhancement of the Josephson current in each 
Andreev levfl. By setting a = L/2 in Eqs. (7), (15), and 
(16) of Ref.El, we obtain 

/± ^ ± 'y^,^ .. Vrcos(^) . (5) 

L + 2e(£;±(0)) ^2' ' 

for symmetrical junction. In short junctions the varia- 
tion of E^{(l)) in Eq. (0) produces an Andreev level cur- 
rent proportional to TH^ while in long junctions Ei^{(j)) ~ 



constant when the barrier transmission T <^ 1. Conse- 
quently, /„ cx Vt in long, low transmission, symmetric 
Josephson junctions, again for the same reasons as the 
resonant enhancement ofjjpsephson current in d-wave su- 
perconducting junctionsEj. Eqs. (^-(j^Lare the same as 
found by Wendin and Shumeiko in Rcf.EZl. 

Fig. |j(b)-(c) shows the Josephson current switching for 
the long, low transmission, symmetric Josephson junc- 
tion. (We fix = 0.3(27r) in Fig. |(b)-(d).) The changes 
in the Josephson current have equal magnitudes until 
the gate voltage approaches the energy gap, similar to 
the ballistic SNS junctionll3. Unfortunately, this switch- 
ing occurs only over a much narrower range of gate 
voltages than in the ballistic SNS junction. Compar- 
ing Figs. @(b)-(c) we see that the range of gate voltages 
where the Josephson current changes from its equilib- 
rium value becomes much narrower as the transmission 
decreases. The VT versus T effect is also clearly visible m 
Fig. ||(c), though none of the 'giant' Josephson currentscZl 
in Fig. H(c) are as large as the ballistic SNS junction in 
Fig. ||(b). The ballistic SNS junction has both the largest 
equilibrium Josephson current and the largest switching 
amplitude of the Josephson current with variation in the 
gate voltage. 

The range of voltages over which the Josephson current 
differs from its equilibrium value in a ballistic SNS junc- 
tion is approximately half of the energy gap, as shown 
in Fig. |j(c). In contrast, the Josephson current in a 
long SIS junction with inversion symmetry differs from 
its equilibrium value only over a very narrow range of 
voltages, as also shown in Fig. ^(c). We can infer from 
Eq. (Q) that the range of voltage over which the 'gi- 
ant' Josephson current of RefEZi occurs is proportional 
to VtA{^o/{L + 2^0)), namely the bandwidth of the 
Andreev level. Fig. 0(d), for the same junction having 
transmission probability T = 12%, shows the differential 
conductance dI]\[/dV along the probe has a peak when- 
ever the probe voltage crosses a new Andreev level. 

B. Asymmetric Junction 



Fig. H shows the (a) Andreev levels, (b)-(c) Josephson 
currents, and (d) gate current in an asymmetric junction 
where a = L/5. The Andreev levels in Fig. |^(a) cor- 
respond to a long SIS junction {L = 6.6^0) where the 
transmission probability T — 100%, 36%, 12% and the 
coupling strength e = 0.1%. The presence of an impu- 
rity removes all of the degeneracy in the Andreev level 
spectrum in an asymmetric SIS junction, though some 
energy levels split more than others. For example, the 
energy gap at = tt and \E\ ~ 0.7 A is much smaller 
than the other gaps. 
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FIG. 5. A long {L > ^o) SIS junction without inversion symmetry, (a) Andreev levels, (b)-(c) Josephson currents, and (d) 
differential conductance along the normal probe. A small Andreev energy gap leads to large nonequilibrium Josephson currents, 
but over a narrow voltage range. Larger Andreev energy gaps produce a nonequilibrium current of nearly the same magnitude 
as the equilibrium current, but over a wider voltage range. Ballistic junctions (T = 100%) again have the largest currents, both 
in and out of equilibrium. 



The small energy gaps present in the Andreev level 
spectrum, even in asymmetric SIS junctions, produce 
Andreev level currents which approach tiie 'giant' Vt 
variation of the Josephson current of Ref.c3. Figs. ||(b)- 
(c) show the Josephson current switching in a long SIS 
junction (L = 6.6^o) having coupling strength e = 0.1%. 
We also fix the phase at = 0.3(27r) in Figs. |5|(b)-(d). 
For the small transmission coefficients in Fig. ||(c) , we see 
that the equilibrium Josephson current decreases much 
more rapidly than the current carried in the Andreev lev- 
els near \E\ ~ 0.7 A (due to the small energy gap near 
\E\ ~ 0.7A). Comparing the asymmetric junction of 
Fig. ||(c) and the symmetric junction of Fig. ^(c) shows 
that the magnitude of the current switching due to occu- 
pation of a new Andreev level can be nearly the same for 
both symmetrical and asymmetrical junctions. There- 
fore, inversion symmetry is not a necessary condition for 
an Andreev level to carry a 'giant' Josephson current. 



The ballistic junction (T ~ 1) again carries the largest 
equilibrium and nonequilibrium Josephson current. 

One further difference between the long SIS junctions 
with inversion symmetry and the asymmetric junctions is 
that the Josephson current switches by different amounts 
in the asymmetric junctions when eV < A, as shown in 
Fig. ^(b)-(c). In an asymmetric junction, magnitude of 
the current flowing through each Andreev level is 
in general different, while for symmetric junctions they 
are nearly the same (as long as \En\ is not too near the 
gap edge). Fig. ||(b)-(c) and Fig. ||(b)-(c) also show that 
Josephson current does not fall exactly to zero when 
probe voltage exceeds the energy gap {eV > A). As 
emphasized in Ref.E3, this is because the leaky Andreev 
levels outside the superconducting gap carry a portion of 
the Josephson current. The normal metal probe there- 
fore provides a means of doing energy spectroscopy of 
the Josephson current, making it possible to measure this 
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'continuum' contribution to the current. Spectroscopy of 
the bound levels is again shown by the differential con- 
ductance diN /dV in Fig. ||(d) for transmission T = 12%. 
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FIG. 6. Andreev levels E (0) and nonequilibrium Josephson current-phase relation /(</>, V) for a junction without inversion 
symmetry (a)-(b) and with inversion symmetry (c)-(d). The current phase relations are similar, except for the larger current 
magnitude in the symmetric junction. 



C. Current-Phase Relation 



The current flow through bound Andreev levels can 
also be obser'ued in the current-phase relation of a Joseph- 
son junctionO In Fig. ^, we show the Andreev lev- 
els and current-phase relation in a long (L — 4^o) SIS 
junctions. The asymmetric junction in (a) and (b) has 
the impurity placed one third of the distance across the 
normal region (a = L/3), while the symmetric junc- 
tion in (c) and (d) has a = L/2. For both junc- 
tions, the transmission probability T = 2.5% and the 
coupHng strength e = 0.1%. We apply gate voltages 
to the SIS junction which intersects an Andreev level, 
namely eV = ±0.414, 0.418, 0.420, 0.425, 0.430A in (a) 
and eV = ±0.465, 0.475, 0.485, 0.495, 0.505A in (c), as 
shown by the horizontal lines in Fig. ^. 

Qualitatively, current-phase relation is similar for both 
types of Josephson junctions. As the quasi-Fermi energy 



sweeps through the a bound Andreev level, the current- 
phase relation changes from sinusoidal, to a half-peijiodic 
relation, and finally to a 7r-phase shifted junctionE3. All 
these changes in the current-phase relatiorustte similar to 
those occurring in ballistic SNS junctions.Bll3, including 
the half-periodic and 7r-phase shifted (sign change of the 
Josephson current) current-phase relations. The current- 
phase relation for the symmetric junction evolves from 
the usual sinusoidal form into a nearly linear variation 
of current with phase more typical of ballistic junctions. 
The half-periodic current-phase relation can be..observed 
as a doubling of the AC Josephson frequencyOc3. 

When we depopulate one of the Andreev levels with 
the probe voltage, we increase the magnitude of the crit- 
ical current. The equilibrium Josephson current is a small 
difference between much larger positive and negative cur- 
rents flowing in equilibrium. Depopulating an Andreev 
level removes some of this current cancellation in both 
the asymmetric and symmetric junction types, increas- 
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ing the Josephson critical current. However, the Andreev 
levels in the SIS junction with inversion symmetry carry 
a slightly larger current (by a factor of ~ 4 in Fig. H). 
There is a small difference between the Josephson cur- 
rent betweeen applying positive and negative gate volt- 
ages not shown in Fig. ||, similar to Fig. g(b) for the short 
junction, again due to a small leakage current from the 
gate. 



IV. CONCLUSIONS 



An additional normal metal terminal weakly coupled 
to a Josephson junction, permits one both to determine 
the Andreev energy levels and to probe the Josephson 
current carried through these bound levels. The differen- 
tial conductance cUn (</>, V) / dV along the normal metal 
probe determines the Andreev level positions and width, 
while changes in the Josephson current I{(f), V) as a func- 
tion of the probe voltage determine the Andreev level 
currents. Setting the probe voltage above the energy gap 
also allows a measurement of the 'continuum' Joseph- 
son current, which flows outside the energy gap. In 
this paper we have studied the Josephson current switch- 
ing and spectroscopy of the Andreev energy levels in a 
three terminal SIS junction, where the normal region of 
the Josephson junction also contains an insulator having 
transmission probability T < 1. The results are qual- 
itatively different from the ballistic Josephson junction 
{T — 1) we considered in an earlier paperli3. 

In a short Josephson junction (L <^ S^q) containing a 
tunnel barrier (T <IC 1), the maximum Josephson current 
switches from the Ambegaokar-Baratoff value {eAT/2h) 
to zero when the probe voltage is approximately the en- 
ergy gap, \eV\ ~ A. In a short junction, there are 
two Andreev levels which carry equal and opposite cur- 
rents. Therefore, populating (or depopulating) both lev- 
els forces the Josephson current to zero. Since the pres- 
ence of a tunnel barrier with T ^ I forces the Andreev 
levels to the gap edge at E ~ ±A, the Josephson cur- 
rent switches to zero when \eV\ ~ A (independent of the 
phase difference). Although the magnitude of the change 
in the Josephson current is larger in ballistic junctions 
(approaching eA/h), the switching voltage in ballistic 
junctions ranges between < eV < A (depending on 

In a long Josephson junction {L ^ ^o) containing a 
tunnel barrier (T <C 1), the details of this current switch- 
ing depend on the scattering potential inside the nor- 
mal region. For the long SIS junction having inversion 
symmetry, a current j^roportional to Vt indeed flows 
through the junctioncZI when one applies certain volt- 
ages to the normal metal probe. (The exact current is 
evp/iL + 2£,{E^))Vt where ^(E^) is the energy depen- 
dent coherence length.) The factor of VT arises because 
the inversion symmetry allows degenerate energy levels 



to exist even in the presence of a tunnel barrier, siifc 
ilar to the midgap states in d-wave superconductors. c3 
Although this nonequilibrium current (proportional to 
Vt) is much larger than the equilibrium Josephson tun- 
neling current (proportional to T), both the changes in 
the nonequilibrium Josephson current and magnitude of 
the equilibrium current are much smaller in SIS junc- 
tions than for ballistic SNS junctions (which approach 
evp/iL + 2^0). 

The bandwidth of the Andreev levels also becomes 
smaller with decreasing barrier transmission in long SIS 
junctions. The range of gate voltages for which one ob- 
tains a nonequilibrium current in the long SIS junction 
is equal to this bandwidth, namely 4\/rA(^o/(i + 2^o))- 
There are approximately {L + 2^o)/27rfo of these en- 
ergy levels, so the total range of voltages over which the 
Josephson current differs significantly from equilibrium is 
AVTA/tt. Reducing the barrier transmission T to max- 
imize the size of the nonequilibrium current simultane- 
ously lowers the range of gate voltages over which one 
can observe this current. In a ballistic SNS junction, 
the Josephson current differs significantly from its equi- 
librium value over a much larger range of voltages A/2, 
namely half of the energy gap. Nonequilibrium effects on 
the Josephson current from the additional normal metal 
probe are much larger and occur over a much broader 
range of gate voltage in ballistic SNS junctions. 

For long SIS junctions which do not possess inver- 
sion symmetry, the energy gaps and currents carried by 
the Andreev levels are in general different, so that the 
Josephson current switches by different amounts when- 
ever the gate voltage populates a new Andreev level. 
The magnitude of the Josephson current switching can 
range between the equilibrium value of the current (pro- 
portional to T) and the larger nonequilibrium currents 
found in symmetric junctions (proportional to Vt). If 
only a very small energy gap occurs near an Andreev level 
crossing, main features of the Josephson current switch- 
ing and spectroscopy of the Andreev levels are qualita- 
tively similar to the symmetric junctions. But because 
the Andreev levels carry different currents in general, the 
cancellation between currents flowing in opposite direc- 
tions from two levels adjacent in energy is almost never 
exact. The Josephson current can differ from its equilib- 
rium value over a much wider range of gate voltages in 
long and asymmetric SIS junctions. 
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APPENDIX A: SCATTERING STATES 



We compute the Josephson current by finding both the scattering state and bound state solutions of the BdG 
equation 



H{x) - n A{x) \ f uix) 

A*{x) -{H*{x)-fi) [v{x) 



E 



u(x) 
v{x) 



For our model Josephson junction we take the Hamiltonian 



H[x) 



2m dx^ 



V{x), 



the scattering potential 

where < a < L, and pair potential 



V{x) = Vs6{x - a) 



A(x) 



Ae^'Pi a; < 
0<x<L 



(Al) 



(A2) 



(A3) 



(A4) 



We first obtain the solutions of Eq. (Al) in a uniform superconductor where A[x) — constant and V{x) = constant. 
When the pair potential and electrostatic potential vary in space, we determine the scattering state solutions of 
Eq. ( Al) by matching quasiparticle wave amplitudes at the potential discontinuities at x = 0, x = a, x and,a: = L 
as shown in Fig. |^. We essentially follow the calculational methods outlined in the Appendices of RefsJa andll^. 



V(x) 



x=a \ \=h J 



N 




FIG. 7. An electron-like quasiparticle is injected from the the left superconductor reservoir. The wave amplitude described 



by -Bo represents the term _B — ^ in Eq. A5. The wave amplitudes in the right superconductor can be obtained interms of A 



and B by matching the boundry conditions at a: = 0, a; = a, a; = 6 and x — L. 

Fig. illusturates an electronlike quasi-particle injected from the left superconducting reservoir. In the Andreev 
approximation, an electron like excitation incident from the left superconductor (x < 0) generates both reflected and 
transmitted quasi-particles described by the wave function 

f uoe**^! \ ^ikeX ^ r ^ „ i£o 1 ( vqc^^^ \ ^ikhx _^ f uqc'^^^ \ ^-ik^x 



valid for (x > < a). In Eq. (|A5|), the first term represents the electronlike quasiparticle injected from the left super- 
conducting reservoir. The term denoted by the coefficient A describes the electron-like particle normally refiected from 
the impurity while the term denoted by the coefficient B describes the Andreev reflected hole-like quasi-particle. The 
transmitted electron-like quasiparticle and the Andreev reflected hole-like quasiparticle in the right superconductor 
(x > L) have the wave function 



\ Vq J \ Uo J 
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The term with the coefficient D describes the hole-like normal reflections from the imp urit y which are first Andreev 
reflected from the right superconductor. The BCS coherence factors uq and vq in Eqs. (A5)-(A6) are 



2ui = l 



\E\ 



and 



2v^ = l- 



\E\ 



The wave vectors fee and kh in Eqs. (|A5|)-(|A6D are 
and 



2m 



(A7) 



(A8) 



(A9) 



(AlO) 



As shown in Fig. |^., quasi-particle waves inside the normal region are grouped into three sets. We obtain the 
quasiparticle wave function in the normal regions (b < x < L) and (0 < x < a) by matching the waves at pair 
potential interfaces x = L and x = respectively. Therefore, for electrons in region (0 < x < a) we have 



+ Ci 



ikhX 







Z /bp oc 





Vn Un I " 



ik 



(All) 



Similarly, in region (b < x < L), matching the wave functions yields 



A. 



Jke{x - b) 



Oj /kh{x~b) ^ 
V 



Jke{x -b) 



Cme-'-~^h{L-b) 



Jkh{x - b) 



For holes in region (0 < x < a) we obtain 
Bi 



1 ^ ^-ik^X ^ jj^ 



1 g-ikhX 



Auqc^'P'^ \ ^-ikeX 





Avo 



-ikfiX 



(A12) 



(A13) 



and in region (b < x < L) we have 



^3 



^-ik,{x -b)^j^ 



) ^-th{x-b) 



DvQe^'t>2/ke{L -b)\ ik^(x _ b) 







~ikh{x - b) 



(A14) 



Since pair potential A is zero in the normal region, the wave vectors k^ and k^ in Eq. (A11)-(A14) are obtained from 
Eq. (^) and ( [A10| ) by setting A = 0. 

The current amplitudes of incoming and outgoing waves at the impurity (x = a) and at the normal-metal side probe 
(x — b) are connected via two separate scattering matrices. The scattering matrix for electrons at x = a is 



rik) i(fce) 



Aie^'^' 

B2 



Similarly, for holes at x = a, the scattering matrix is 



C,e^ha\ ^ ft*{h) r*Ckh) 



D2 



r*{kh) t*{kh) J \Die-'^h 



-ikhCL 



(A15) 



(A16) 
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The current transmission and reflection amplitudes for the point scatterer matrices in Eqs. ( Al5)-( A16| ) are 

1 



and 



r{k) 



1 + i{mV^/n^k) 

-iimVs/n^k) 
l + i{mVjh^k) 



(A17) 



(A18) 



At X = 6 the quasiparticle wave amphtudes are coupled by another scattering matrix to the normal- metal side probe. 
Since we assumed that the normal-metal probe only breaks the phase of a quasi-particle, but not its momentum, we 
conceptually view the normal probe as two seperate leads. We couple right-moving quasi-particles only to lead 3 while 
the left-moving quasi-particles to lead 4. Therefore, the scattering matrix becomes 



/ 

VT^e 
^^e 

V -VT~^ 
















(A19) 



The scattering matrix S from Eq. ( A19 ) relates the current amplitudes for electrons in normal region (a < x < b) 
and normal region (b < x < L) as 



/ ^3 \ 

\ Bl3 J 



(S) 



/ ,^3 

A^^ikeib-a) 


V 



(A20) 



The current amplitudes for holes, in the same regions are related as 



^3 

V C^L3 J 



(S) 



V 



^3 






(A21) 



By combining Eqs. (|A5D, (|A^), ( |All|) , (|A12[ ), ( |A13|) , (|A14| ), ( |A15|) , (|A16D , (|A20|), and ( |A2lD we obtain matrix 
equations for the coefficients C and D as 



tike) r{ke) \ ( e^(ke-kh)c 



(1 _ ^)e'^(t>e-i{ke - kh){a - L) 



Wo 



uoj \r*ikh) t*{kh) 



X 







(A22)\ 



J{kf, - kh)a 



(uo_ Vo \ pikhtt 




1 Q\ C 



1 



A' 



(A23) 



and 



t*{kh) r*{kh) 





t{ke) r{ke) \ ( e^(ke - kh)a 

r{ke) tike) 





J{ke - hi)a 



X 







[I _ ^y~i{ke - kh){a - L) 



B' 
D' 



(Uo _ Vq_ 



Jk/iU 



1 0\ f B' 



1 



D' 



(A24) 



In Eqs. (A23), the C", A\ B' , D' are related to C, A, B, D through 
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(A25) 



A! = Ae 



-ikpa 



(A26) 



and 



D 



B' = Be^^h"' , 
jj^ikh{L - a) 



(A27) 



(A28) 



Finally, by inverting the matrix equations (A23)-(A28), we can directly obtain the scattering amplitudes as 



1 _ ( ^ I (1 _ e)e*'^e*(^- " '"^^^ VT~e 1 - 



F(0,e,T) 



(A29) 



and 



Dl 



""0/ y \uo. 



The resonant denominator F{(f>, e,T), whose zeroes give the position and width of the Andreev resonnaces, is 



(A30) 



uo 



F(0,e,T) = l-2T(l-e) ( — ) e*(^'= " cos - i? | — | e 



^2i{ke-kh)a 



(A31) 



For any other quasiparticle, injected from either superconducting contacts or "conceptual" normal probes, one can 
obtain the wave-function amplitudes by calculations similar to the ones described above. In Appendix B we outline 
the derivation of the electrical current flowing into the right superconductor from those wave-function amplitudes. 



APPENDIX B: ELECTRICAL CURRENT 

The electrical current flow in the right superconductor due to the injection of a quasiparticle from any lead i 

I2 = J2 ('^" + J^ f,^2 /f - E ■ (Bl) 

fe;g/3 fc:ij/3 

We use the indices q — 1,2,3,4 as the lead numbers, k as the wave number of the injected quasi-particle, and 
/3 = (e or h) denotes the injection of electron-like or hole-like quasi-particles. The J„ and J„ are the Schrodinger 
currents carried by the waves u and u, namely, 

J„ = {eh/m)lm{u*{x)Vu{x)} , (B2) 

and 

= ieh/m)lm{v*{x)Vv{x)}. (B3) 
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The last term in Eq. (Bl), ^ ^..^^ {Jv)q^2 ~ is called 'vacuum current' and is zero in DC problems. 

We can now define the transmission coefficients r^_^2i which give the ratio of electrical current out to particle 
current in, as 



where vp = hkp/m. Equation (Bl) then reduces to 



k-qP 



which is the standard Landauer-Biittiker form. The Fermi factor in Eq. (B5) is 



(B4) 



(B5) 



(B6) 



where eV^ is the effective biasing voltage applied to the qth lead. Since we ground both superconducting leads, we 
have eVf = = 0. The effective biasing voltages applied to the normal-metal leads are 



eVi = eVi = eV, 
eV^ = elf = -eV 



(B7) 
(B8) 



We obtain the current flow into the right superconductor by converting the sum over k in Eq. (B5) into an integral 
over the injected energies, namely 

-oo J A 



vfN+{E)[T^ 



-n^2]f{E~eV)dE 



/OO 
vfN+{E)[t!^^^ + T':^2]f{E + eV)dE , 
-oo 



where the superconducting density of states N^{E) related to the normal density of states N^{E) by 

1 



Nt{E) 



K-vl\ 



K{E) . 



(B9) 



(BIO) 



Combining the identity vfN+{E) ~ l/h, along with Eqs. (B9)-(B10), we obtain 

f-A 



h 



,e,V,T)^-{ 



+ TL2 + + T^\2]fiE)dE 
I [Tl^2 + TU2]fiE - eV)dE 

/OO 
[T^^2+n^2]f{E + eV)dE , 



(Bll) 



The transmission c oeffi ci ents T in Eq. (Bll) are related simply to the scattering amplitudes calculated in Appendix 
A. Using Eq. (^), (^, (^ and (^ we obtain Tf^2^ for example, as 



\Dl_ 



(B12) 



The other transmission coefficients and wave-function amplitudes can be obtained by calculations similar to the one 
in Appendix A. Hence for a hole-like quasiparticle injected from the left superconductor we have 



T'l^2 = \Ci^2?-\Dt 



■2I ' 



(B13) 



where 
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F(,.e.r) • C"*) 

and 

= '-^ • 

For an electron-like quasiparticle injected from the right superconductor we find the transmission probability 

Ti^2 = -1 + |C|^2l' - 1^1^21' - (4«o^^o - 2^)Re{D'^^,} - ^ + 4«o , (B16) 

where 

2 - -\/ / \ 2\ 



and 



^2^2 = 7 ^ — —■ (B18) 

For a hole-like quasiparticle incident from right superconductor, we find 

T2t2 = 1 - \C^^2\' + 1^2-^2 1' + (4wot'o - 2^^)Re{D^^^} + ^ - (B19) 



where 



= ^ — ^-^ ^ , (B20) 

F{(j), €,T) 



and 



il-(^) e2^(^--^'')« l-(^) (l-e)e2^^''(^-«) 
D^^2 = r* ^ • (B21) 

For an electron-like quasiparticle injected from the normal probe in channel 3 we have 

^3-^2 = \Ct^2? - 1-^3^21^ 1 (B22) 

where 

h-f^V e^^{ke-kh)a \ ^/J^^^^^^rj^^/^^T^^-iheML + b - 2a) 
Cl^2 = r- — =7 , (B23) 

and 
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yR{l - e)e*'^^e - kh){L - a) ^ 2-g#g«(^e - kh)a _ _ ^^giike - kh){L + a) 

uoF{<p,e,T) 

For a hole-like quasiparticle incident from the normal probe in channel 3 we have 

r3t2 = |^3^2p-|C3^2p, 

where 

^3^2 = ( -) V|«oP + \vo\'V7{l^)e-'~^h{b - a)/ke{L - a) 



Uo 



uoF{(l), e,T) 

and 

2 



1_ 1^^^ e2iCke-kh)a\ ^|„q|2 + |^Q|2y^(r3^e-ifcft(i. Hr 6 - 2a) 



n^i = r* 

uoi^(</.,e,T) 

For an electron-like quasiparticle injected from the normal probe in channel 4 we have 
where 

C|^2 = V\uo\^ + \vo?V~e 
^1 _ [i?e^*(^e ~ ^'')" + T(l - e)e*'^i?e*(^e ~ ^/i)-^-] j ^-i(f)2 ^ike{L - h) 

uoF{<l>,e,T) ' 

and 

DU2 = r* (-) Vl«oP + - e) 

^1 _ g2i(^e - kh)a^ ^-i(j)2^ike{L - b)^-2ikh{L - a) 

woi^(<^,e,T) ' 
For a hole-like quasiparticle incident from the normal probe in channel 4 we have 

where 



CL2 = M ^ ) Vl^^oP + KlW - e) 



1 _ ( g2i(A;e - kh)a^ g-ikh{L - b)^2ike{L - a) 
uoF{(t>,e,T) ■ 
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/I 



icj) i{ke - k,i)L-i I p-ikh{L - h) 



uoF{(j),e,T) 



(B33) 
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